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We present the results of simulations of the ultra-high energy cosmic ray (UHECR) propagation in the Galactic 
magnetic field (GMF). Different assumptions on the large-scale GMF structure and/or primary particle lead to 
distinctly different deflection patterns of the highest energy cosmic rays (CR). The GMF modifies the exposure 
of an UHECR experiment to the extragalactic sky. We estimated these effects for the Pierre Auger experiment. 
Further forward-tracking studies under plausible UHECR sources scenarios will allow for direct comparison with 
the observed correlation between the nearby active galactic nuclei (AGN) and the highest energy Auger events. 



1. Introduction 

Magnetic fields in the Milky Way are inves- 
tigated via Faraday rotation of the polarized 
light from pulsars and extragalactic (EG) sources, 
and through measurements of the Galactic syn- 
chrotron radiation [12 . The observations reveal 
the presence of a large-scale (LS) Galactic mag- 
netic field (GMF), that in the first approximation 
follows the spiral arm structure. The LS field ex- 
tends above and below the Galactic disk on a kilo- 
parsec scale and forms a kind of halo. A random 
GMF, of the strength similar to that of the reg- 
ular component, but with the coherence length 
of only ~ 50 pc, has also been observed [1 . The 
average strength of the total magnetic field near 
the Sun is about 6 jaG. However, details of the 
field distribution are poorly known ^3p4p5j . 

The Pierre Auger Observatory |6] provides a 
new and independent way of studying cosmic 
magnetic fields, by collecting CR events above 
10 EeV = 10^^ eV with unprecedented statistics 
and data quality. Recently, Auger observed a 
significant correlation (over angular scales < 6°) 
of the arrival directions of CR above 60 EeV 
with the locations of nearby AGN [7 8 , and a 
strong steepening of the CR flux above 40 EeV [9 . 
Both observations are consistent with the stan- 
dard UHECR astrophysical acceleration scenar- 
ios and thus represent an important step towards 
the "charged particle astronomy" . 



The elaboration of relevant analysis methods 
for the UHECR astronomy requires detailed in- 
vestigation of the cosmic ray propagation in the 
GMF. We present in this paper results of such 
studies, performed in the light of the observed 
AGN correlation. We used the standard method 
of CR backtracking. A large number of CR events 
has been simulated using parameters described 
below, and propagated under three distinctive 
large-scale GMF models. We present the result- 
ing magnetic deflections and modification of the 
extragalactic exposure due to the LS GMF. 

2. Implementation of CR backtracking 

To obtain the UHECR trajectories, we imple- 
mented the integration of the equations of mo- 
tion of an ultrarelativistic particle in the magnetic 
field using the Runge-Kutta 5*^ order scheme 
with the adaptive step size control [10 . The 
particles have been followed till the galactocen- 
tric distance of 20kpc (the Galaxy "border"), be- 
yond which the GMF strength is supposed to be 
negligibly small. An accuracy level (the accepted 
truncation error) of 10~^ has been adopted, for 
which the changes in the backtracked directions 
become negligible with respect to the chosen step 
in angular separation (imax from the selected as- 



trophysical objects (see the section 6.1). To avoid 
the "numerical dissipation" of the CR energy [11 , 
we preserved the absolute value of the CR veloc- 
ity vector during propagation. 
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(a) Bisymmetric even parity spiral field (b) Axisymmetric odd parity spiral field 



Figure 1. Cosmic ray deflection patterns for protons under Models I and 11. Arcs of great circles join 
the simulated arrival directions on Earth, and the backtracked ones at the Galaxy border (arrow heads). 
Every 1000*^ event is shown. The dashed line denotes the supergalactic plane (SGP). The solid lines 
show the equal integrated Auger SD exposure sky regions within the detector field of view for Oz < 60°. 



3. Parameters of simulated CR events 

Energies of the simulated events were 
bounded between 40EeV and 150 EeV, and dis- 
tributed according to the Auger measurement [9j , 
as a power law 

Two options for the angular distribution of 
the events have been considered: 1) Uniform over 
the whole sky (10^ events), and 2) According to 
the exposure of the Auger Surface Detector (SD) 
for zenith angles 6^ < 60° [12 13 (10^ events). 

Four pure mass compositions have been con- 
sidered: protons, carbon, silicon, and iron nuclei. 

4. Considered large-scale GMF models 

We have chosen amongst many available large- 
scale GMF models three typical ones with distinct 
qualitative differences. Two models are the spi- 
ral disk field models of bisymmetric even parity 
(hereafter Model I), and axisymmetric odd par- 
ity (Model II) from the paper [IT by Harari, 
Mollerach and Roulet (HMR). 

Model III is a modification of the model [15] 
by Prouza and Smida (PS), that has been pro- 
posed by KachelrieB et al. In addition to 
the bisymmetric even parity spiral field, of the 
structure similar to the one in Model I, it fea- 
tures two additional large-scale halo GMF com- 
ponents: toroidal azimuthal field above and below 



the Galactic disk, and poloidal (dipole) field. 

5. Magnetic deflections of cosmic rays 

The cosmic ray deflection patterns obtained us- 
ing the three selected GMF models are very dif- 
ferent, and have trends, typical for each model. 
Fig. [T] shows deflections of primary protons for 
Models I and II. In the bisymmetric even parity 
field, CR arrival directions at the Galactic bor- 
der point in general to higher Galactic latitudes 
than observed on Earth. In the odd parity field, 
the backtracked arrival directions in each Galac- 
tic hemisphere are shifted to the poles, and no 
cosmic rays are coming from the directions close 
to the Galactic plane. The presence of dipole and 
toroidal fields in Model III makes the deflection 
pattern even more complex [17]. 

The deflection values for the three GMF mod- 
els, and the Auger SD exposure are summarized 
in the table [T] by means of percentiles at 50% 
(median), and 95.45% of the c.d.f. The median 
deflection values scale well with the atomic num- 
ber Z of primary nuclei. They are larger for the 
HMR models than in Model III, due to the more 
important halo field extension. 

The simulated events allow to roughly estimate 
the corresponding deflections expected from tur- 
bulent magnetic fields [17 . Compared to the de- 
flections from the large-scale field, the former ones 
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have to be much smaher, which is confirmed by 
simulations including turbulent field [18^, and can 
be neglected at the first approach. 

6. GMF effects on the EG exposure 

(De-) focusing effects of the Galactic magnetic 
field make the correspondence between the CR 
arrival directions on Earth, and the EG sources 
contributing to the CR fiux entering the Galaxy 
non-trivial jl4|16|19 . We have constructed maps 
of parameter A = log^g of ratio of the number 
of CR arrivals on the halo border to the one on 
Earth, using the HEALPix equal area celestial 
sphere pixelization [20 . These maps are shown 
on Fig. [2] for the case of primary carbon nuclei. 

Due to these effects, some regions of the sky 
have increased probability to contribute to the 
cosmic rays observed on Earth, and some others 
are disfavored. The non-uniform character of the 
mapping strengthens in the case of heavy compo- 
sition [17 , so that the regions that would effec- 
tively contribute to the CR fiux on Earth repre- 
sent only a small fraction of the Att solid angle. 




6.1. Correlation scan using backtracked di- 
rections 

To quantify the lensing GMF effects on the EG 
exposure for the Pierre Auger Observatory, we 
performed a correlation analysis using backtracked 
arrival directions of simulated events, and 694 
AGN at redshift z^ax < 0.024 from the VCV cat- 
alogue [21 . The simulated CR have been divided 
sequentially into samples with the same number 
of events (81) above 40EeV as in the Auger data. 
The employed set and ranges of parameters were 
also identical to the ones from [7fS^ . 

We will focus here on the minimum probabil- 
ity values, obtained during the scan. The c.d.f. 
of decimal logarithm of the corresponding cumu- 
lative binomial probability Pmin of reaching this 
level of correlation under isotropy for Model I 
are shown on Fig. [3j The level of the minimum 
probability reached in the Auger data is also in- 
dicated. The results of the correlation scan for 
the assumed LS GMF models and primary mass 
compositions are summarized in Table [T] 

Since for the different GMF model/primary 



Figure 2. Sky maps of the A values for Models I 
(upper), II, and III (lower), and carbon nuclei. 
Earth UHECR detectors are (almost) blind to the 
regions in gray. The dashed line denotes the SGP. 



mass assumptions the backtracked directions cor- 
relate with the catalogue objects in particu- 
lar privileged regions in the sky, our scan re- 
sults depend strongly on those assumptions. For 
Model III, the scanned probability minimum is 
significantly less deep than for the two other 
spiral-field-only models for any assumed primary 
mass, except for protons, where one obtains 
nearly the same level of correlation for all mod- 
els. This model is clearly less compatible with the 
observed correlation with the nearby AGN [7 8 , 
unless the UHECR fiux contribution from these 
objects (or objects with similar spatial distribu- 
tion) is highly non-uniform. 
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loglO(Pmin) 



Figure 3. Scan results (c.d.f. of log^^g ^min) of 
the simulated CR samples, backtracked under 
Model I and 4 masses. The dashed line shows 
the value from the Auger scan [7 8 on Earth. 
C.d.f. for protons (shown by symbol at the me- 
dian value) does not intersect this line. 
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HMR axisymmetric odd parity model 
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-9.6 
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version by KachelrieB et al. 


P 


3.0 


30.1 


-4.2 
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-1.5 
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17.6 
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-4.3 


-2.5 


-1.4 
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92.0 


-4.9 


-3.0 


-1.7 


Fe 


58.3 


124.7 


-4.2 


-2.6 


-1.5 



Table 1 

Magnetic deflections (in °), and log^^g ^min values 
at the indicated percentiles of the c.d.f., for the 
assumed GMF models and primary mass. 



7. Conclusions and outlook 

Our studies of the UHECR propagation in the 
LS GMF show that the UHECR picture observed 
on Earth is sensitive to the assumptions on the 
field structure and/or primary CR composition. 

The exact defiection value in the regular com- 
ponent of the field depends strongly on the arrival 
direction on Earth of a cosmic ray, with the cor- 
responding position angle of defiection differing 
from one assumed field distribution to another. 
To discriminate between GMF models, additional 
hints can be provided by the aligned structures of 
events coming from a UHECR source [22 . 

Though the reconstruction of the field is eas- 
ier in the case of light primary mass composition, 
in the case of heavy nuclei the lensing effects of 
the Galactic field on the exposure bring stronger 
constraints on the list of potential UHECR source 
candidates. The presented analysis of the corre- 
spondence between the arrival direction distribu- 
tions on Earth and at the Galactic border will be 
complemented by the forward-tracking of cosmic 
rays for a number of plausible UHECR sources 
scenarios. This will allow for direct comparison 
with the observed AGN correlation. 
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